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The transverse mass spectra and mid-rapidity yields for Hs and fis plus their anti-particles are 
presented. The 10% most central collision yields suggest that the amount of multi-strange particles 
produced per produced charged hadron increases from SPS to RHIC energies. A hydrodynamically 
inspired model fit to the spectra, which assumes a thermalized source, seems to indicate that these 
multi-strange particles experience a significant transverse flow effect, but are emitted when the 
system is hotter and the flow is smaller than values obtained from a combined fit to 7T, K, p and As. 

PACS numbers: 



In heavy ion collisions we aim to investigate nuclear 
matter under extreme conditions of pressure and tem- 
perature which is expected to lead to the creation of 
deconfined partonic matter, the Quark Gluon Plasma 
(QGP) lj. The study of strange particles production 
is thought to yield information on the collision dynamics 



from the early stage to chemical and thermal freeze-out. 
The production of strangeness through partonic interac- 
tions, mainly gg — > ss, is expected to dominate over that 
by hadronic scatterings 0, HI- Since at RHIC energies 
the initial gluon density is expected to be much higher 
than in lower-energy heavy ion collisions 0, |]| , strange 
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particles, especially multi-strange particles are expected 
to be more sensitive to the partonic dynamics of the early 
stage of the collision. 

Chemical freeze-out is defined by the temperature, T c h, 
at which inelastic collisions cease and the relative parti- 
cle ratios become fixed. With no strange quarks in the 
initial state it is possible for the system to freeze-out 
chemically before absolute strangeness chemical equilib- 
rium has been reached. Statistical models have had much 
success in reproducing the particle ratios both at lower 
energies (e.g. HH||) and at RHIC (e.g. HG3). Some 01 
these models, which assume chemical equilibrium, allow 
for the possibility of non-equilibrium in the strangeness 
channel by the introduction of a strangeness saturation 
factor, 7 S 0>E3- This factor reaches unity when com- 
plete strangeness saturation, or absolute chemical equi- 
librium, is present in the system. A common T c h of 170 
MeV |2| is indicated by fits to measurements at 130 GeV 
of anti-baryon to baryon ratios as well as particle ratios 
involving only mesons and non-strange baryons. This 
temperature is remarkably close to that of the phase tran- 
sition calculated with lattice QCD Re-calculating 
the fits with the inclusion of the multi-strange baryons 
allows a more sensitive test of statistical models to be 
made. 

A comprehensive analysis of the transverse momentum 
or transverse mass spectra allows us to probe the possi- 
bility of thermal equilibrium of the system. The high 
mass of the multi-strange baryons means they are more 
sensitive to the size of the transverse flow of the system. 
It has previously been suggested that thermal freeze-out 
occurs much earlier for Q and S than for lower mass par- 
ticles due to the predicted low scattering cross section 
of fl and S I n the limit of vanishing cross sections 
this would mean that these particles are emitted almost 
directly from the phase boundary of the hadronizing fire- 
ball. However, these cross sections have not actually been 
measured and may not be negligible. If they were signifi- 
cant, combined with the now sizeable particle production 
at these energies, the fl + X elastic collision rate would 
delay the Q, kinetic freeze-out, resulting in these particles 
carrying the same flow velocity as the pions and kaons. 
By comparing different thermal freeze-out scenarios to 
the data (e.g. G3E1E3) it may be possible to deter- 
mine the rate of build-up of radial flow. 

The data presented here were taken with the STAR 
detector in Au+Au collisions at y/s NN = 130 GeV. The 
main tracking detector is a large cylindrical Time Pro- 
jection Chamber. A Central Trigger Barrel measuring 
the produced charged particle multiplicity around mid- 
rapidity plus two Zero Degree Calorimeters essentially 
measuring neutral spectator energy were used for trig- 
gering lid . In this paper the centrality of the collisions 
is determined from the measured mid-rapidity negative 
particle multiplicity. The data were divided into three 
centrality classes corresponding to 0-10%, 10-25% and 
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FIG. 1: m± spectra for (a) H - , (b) E + and (c) Q as a function 
of centrality. Scale factors have been applied to the spectra 
for clarity. Points are drawn at the bin center. The horizontal 
bars indicate the bin size. The dashed curves are Boltzmann 
fits to the spectra. The solid curves are hydrodynamically 
inspired model fits to the most central H~and H + spectra. 
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0.27±0.02 
0.22±0.02 
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TABLE I: 



H + and O, 



Fit parameters for the m± spectra of the 3 , 
The data represent the 0-10%, 10-25% and 25- 
75% centrality bins with h~=dN h ~ /drj\ |^|<o.5- Only statisti- 
cal and p± dependent systematic uncertainties are presented. 
The p± independent systematic uncertanties are 10%. 



25-75% of the total hadronic cross section as described 
in [l7j. Multi-strange particles are identified via their 
decay modes 3 — > A + ir and £1 — > A + K with the 
subsequent decay of A — > p + n. The tertiary A ver- 
tex is identified by selecting positive and negative tracks 
that are consistent with an origin at the decay of a hy- 
peron some distance from the primary collision point [18j . 
The secondary vertex of the decay is located in a simi- 
lar fashion by combining the previously identified A with 
a charged particle. Simple cuts on geometry, kinemat- 
ics and particle identification, via specific ionization, are 
applied at each step to reduce the background due to 
the high multiplicity ^3 ■ The momenta of the daughter 
particles at the decay vertex are then combined to cal- 
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FIG. 2: (a) The kinetic freeze-out temperature vs transverse flow velocity for the hydrodynamically inspired model fits to the 
m± spectra. The 1, 2 and 3 sigma contours are shown. Solid curves are for a simultaneous fit to the H~and E . Dashed curves 
are a separate fit to the STAR n, K, p and A data. The diamonds represent the best fit in both cases, (b) Mean transverse 
momenta for identified particles vs particle mass (see text for details). The band results from the three sigma contour of the 
hydrodynamically inspired model fit to the ir, K, p and A data and the dashed curve is for Tf = 170 MeV, < f3± >= 0. 



culate the parent particle kinematics. The peaks in the 
invariant mass plots have an average signal to noise ratio 
of 0.74, 0.78 and 0.86 for the S~, S + and f2 respectively 
in the 0-10% centrality bin. The signal to noise ratios 
are calculated for ±15 MeV/c 2 about the expected mass. 
The statistics of the fi~and fl signal are not sufficient 
to allow a separate measurement of the spectra of each 
particle. Hence f2 refers to f2~plus fl + . The momentum 
integrated f2 + /fi~ratio for the top 11% most central data 
is 0.95 ± 0.15(stat) ± 0.05(sys) 0. 

The invariant mass distributions are histogramed in 
transverse mass, m± = \J°p\ + to 2 , and the signal ex- 
tracted for each bin. The raw yields are calculated from 
the invariant mass distributions by counting the entries 
within ±15 MeV/c 2 about the expected mass and then 
subtracting the background. The background is esti- 
mated by sampling two regions on either side of the peak. 
The raw yield in each m± bin is corrected for detector 
acceptance and reconstruction efficiency by the Monte 
Carlo technique, where simulated particles were embed- 
ded into real events. The data cover | y |< 0.75, where 
efficiency and acceptance studies have shown the correc- 
tions to be constant. The total correction factors for the 
S(f2) are 0.2%(0.04%) for the lowest m± bin rising to 
4.0%(0.5%) for the highest bin. 

Fig- UJ shows the invariant m± spectra as functions of 
centrality for the S~and S + , and the for the 10% most 
central data. A portion of the systematic uncertainties 
are a function of p± , therefore this uncertainty is added 
in quadrature, on a bin-by-bin basis, to the statistical 
one, yielding the vertical bars in Fig.^ We estimate the 
remaining systematic uncertainties to be 10% on both 
the extracted invariant yields and slope parameters, the 
major source of which lies in the misrepresentation of the 
embedded Monte-Carlo to the data resulting in a system- 



atic uncertainty in the efficiency calculation. These were 
obtained by exploring the dependence of the invariant 
yields and slope parameters to changes in the cuts phase 
space (more details can be found in [2£j, 0, H3)- The 
weak decay feed-down of fl on S is estimated to be less 
than 2% and thus is neglected. 

Table |I] shows the results of fitting the m± spec- 
tra by an exponential (AEe~^ m± ~ m ^ TE ) and Boltzmann 
(ABm±e~( m± ~' m ^ TB ) functions. Both functions give a 
good representation of the data but the Boltzmann fit 
gives a slightly better y^/dof. The inverse slopes of 
the S particles are the same within uncertainties and 
show no apparent increase over the measured centrali- 
ties. The yields per unit of rapidity are extracted by 
integrating each function over the entire m± range. The 
measured S spectra correspond to ~75% of the total yield 
and the fl to ^66%. The S~and S yields as function 
of ^i~=<A/V; l -/d7/||, ) | < o.5 appear linear; such behavior was 
reported for the A hyperon at RHIC |l8(. 

The 10% most central m± spectra are also fit by a hy- 
drodynamically inspired function |23| . In this model, 
all considered particles are emitted from a thermal ex- 
panding source with a transverse flow velocity < 8\ > 
at the thermal freeze-out temperature Tf D . As in |24| a 
flow velocity profile of f3±(r) = /3 s (r/i?) 5 is used, where 
R is the maximum emission radius. The solid lines of 
Fig. |2Ia) show the one, two and three sigma contours 
for Tf vs < /3j_ > for the fit to the S~and S + data 
combined, with the diamond indicating the best fit so- 
lution (T fo = 182 ± 29 MeV, < X >= 0.42 ± 0.06 c, 
X 2 /dof — 13/15). Also shown, as the dashed lines, are 
the one, two and three sigma contours for a combined fit 
to the STAR 7T, K, p, and A data |H El , and 
the marker is the optimal fit location. Clearly the results 
for the two data sets do not overlap indicating that the S 
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baryons, within this approach, show a different thermal 
freeze-out behavior than n, K, p, and A. The current f2 
statistics do not allow to distinguish between an early de- 
coupling or a common freeze-out with the lighter species. 
Fig. |2fb) shows the mean p± for these particles calcu- 
lated from the functions which best reproduce each m± 
spectrum (Bose-Einstein for n, exponential for K, hydro- 
dynamically inspired function for p and Boltzmann for A, 
S and f2). The error bars on the experimental points are 
statistical and systematic uncertainties added in quadra- 
ture. The band represents the model prediction based on 
the three sigma contour for the fit to the STAR ir, K, p, 
and A data while the lower dashed curve shows the pre- 
diction for Tf o =170 MeV, <(3±>~0, i.e. a system where 
thermal and chemical freeze-out coincide and no trans- 
verse flow is developed. From Fig. EJa) it is likely that 
the S baryons prefer a hotter thermal freeze-out temper- 
ature parameter when compared to that resulting from 
fits to the lower mass particles. This is consistent with 
the calculated < p±> of the 2 and Sis being below the 
solid band in Fig.EJb). 

If, as indicated by this fit, the 5 thermal freeze-out 
occurs in conjunction with that of chemical freeze-out, 
T c h ^170 MeV, this could be an indication that a sig- 
nificant fraction of the collective transverse flow has al- 
ready developed at/before chemical freeze-out, probably 
at a partonic stage. Two alternative models have also 
been proposed for describing the RHIC spectra. The 
first [lj| assumes that all transverse radial flow is devel- 
oped at/before T c h=Tf Q =165 MeV. The apparent soften- 
ing of the lighter mass spectra is then due to contamina- 
tion from resonance decay products. A hydrodynamical 
approach [27J is used in the second scenario where the 
particle mean free paths are assumed small until ther- 
mal freeze-out around Tf o =110 MeV. These large cross- 
sections result in even the £7 developing a sizeable fraction 
of its radial flow after hadronization. 

The S + /ft~ and A/ft, - ratios for the most central data, 
as shown in Fig. EJa), increase from SPS energies |2^| to 
RHIC, whereas the S~ jh~ ratio stays constant and the 
A/ft - decreases. When discussing the baryon ratios the 
interplay of increased strangeness production and reduc- 
tion in the net-baryon number, which was significant at 
the SPS, has to be considered. The proximity of the net- 
baryon number to zero at RHIC is reflected in the fact 
that the ratio of S~/ft~ is close to S + /ft~. The reduc- 
tion in the net-baryon number has a larger effect on the A 
than on the E - , as seen in Fig^a), and thus creates the 
observed rise in the S~/A ratio (Fig.|3fb)). It is interest- 
ing to note that the S + /A ratio is a constant from SPS 
to RHIC indicating that the scale of the multi-strange 
enhancement is the same for singly and doubly strange 
baryons. 

A fit to the reported 0-10% centrality particle ratios 
from STAR B El |H H i3 , including the multi- 
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FIG. 3: (a) ~ET ', E + , A and A to ft~ ratios and (b) H~/A 
and H + /A ratios for the most central data as a function of 
y/s NN . The solid lines indicate the statistical uncertainty 
while the caps indicate the statistical and systematic uncer- 
tainties added in quadrature. Some ratios are slightly shifted 
along the x-axis for clarity. 



Ratio (xlO -3 ) 


Data 0-10% 


Stat. 


Non-Equil. HIJING/ BB 


E-/h- 


7.9±0.4±1.0 


7.7 


7.6 


5.1 


E + /h- 


6.6±0.3±0.8 


6.5 


6.1 


3.0 


n /h~ 


2.2±0.5±0.4 


2.9 


2.8 


0.29 


E-/A 


193±18±27 


148 


190 


171 


E + /A 


219±21±31 


163 


207 


142 



TABLE II: Ratios from the 10% most central data compared 
to predictions from a statistical model[10j a non-equilibrium 
model[3(j and HIJING/BB 31]. Uncertainties shown are first 
statistical and then systematic. 



strange particle measurements, using the thermal model 
described in results in a /dof of 8.5/9, a T c h 
of 181±8 MeV, light quark and strange quark chemical 
potentials of 11.7±0.6 MeV and 0.9±1.6 MeV respec- 
tively, and 7 S =0.96±0.06. The fact that 7 S is equal to 
unity, within errors, when used as a free parameter in 
the model indicates that a saturation of strangeness pro- 
duction has been achieved in the most central collisions. 
Table [H] shows particle ratios for the 0-10% centrality 
bin compared to values from three different models: a 
fit corresponding to the statistical model used above, an- 
other corresponding to a statistical model that allows 
for chemical non-equilibrium via an over-saturation of 
both the light quark and strangeness phase spaces |30| 
and the prediction of the event generator HIJING/Bi? 
vl.O which uses a gluon-j unction mechanism to enhance 
the transport of baryons to mid-rapidity [3l| . All ratios 
are well reproduced by both statistical models, indicat- 
ing that the strangeness production is duplicated well by 
such approaches. HIJING/B-B, however, fails to predict 
the multi-strange to ft~ ratios, in particular that of the 
f2, which it under-predicts by nearly an order of magni- 
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tude. It has previously been shown that HIJING suc- 
cessfully predicts the mid-rapidity total charged particle 
yields [32} suggesting that the entropy is reasonably well 
reproduced by this model. The addition of the gluon- 
junction mechanism, which was necessary to replicate the 
small net baryon yields at RHIC |3^| , does not sufficiently 
enhance the multi-strange baryon yields suggesting that 
a different physics mechanism is necessary to model the 
strangeness production. At SPS energies the introduc- 
tion of final state interactions helped to account for the 
observed hyperon enhancement but failed to reproduce 
the overall strangeness production (K + /ir + ). The intro- 
duction of strong partonic interactions in the initial state 
was needed to account for both the hyperon and overall 
strangeness production at the SPS |3J|. 

In summary, the total yield of multi-strange baryons 
per h~ in Au-Au collisions is increased compared to that 
at the top SPS energies. A chemical analysis of the data 
indicates that for central collisions the strangeness phase 
space is now saturated. A fit of the m± spectra to a 
hydrodynamically inspired model suggests the S baryons 
thermally freeze out of the rapidly expanding collision 
region at a hotter temperature, which is close to that of 
chemical freeze-out, and with a smaller transverse flow 
than the lighter particle species. This suggests that they 
decouple at an earlier stage of the collision and thus probe 
a different dynamical region, but one at which a sizeable 
fraction of the transverse flow has already developed. All 
these observations are compatible with the early stage of 
the collision being driven by partonic interactions. 
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